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Abstract. Precise radial velocity observations of the 
nearby M4 dwarf G1 876 with the Observatoire de Haute 
Provence 1.93 m telescope and the new 1.20 m swiss tele¬ 
scope at la Silla indicate the presence of a jovian mass 
companion to this star. The orbital fit to the data gives 
a period of 60.8 days, a velocity amplitude of 246m.s“^ 
and an eccentricity of 0.34. Assuming that G1 876 has a 
mass of 0.3 Mq, the mass function implies a mass for the 
companion of 2/sin i Jupiter masses. 
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1. Introduction 

The still recent discovery of the first extrasolar planet, 
around 51 Peg (Mayor & Queloz, 1995), has since been 
followed by many more. The count currently runs to 11 
very low mass companions (Marcy & Butler, 1998; Queloz, 
1999), with minimum masses (M sini) which range be¬ 
tween 0.5 and 10 times the mass of Jupiter. Asides from 
their Jupiter-like masses, which largely reflect the sensi¬ 
tivity threshold of current radial velocity programs, the 
known extra-solar planets are a very diverse class. Some of 
them have large eccentricities when others have nearly cir¬ 
cular orbits, and their periods range between 3.3 days and 
4.4 years. Giant planets can thus have very much shorter 
periods than in our solar system, which clearly does not 
represent the only possible outcome of planetary system 
formation and evolution. 

To date on the other hand, planets have mostly been 
looked for around solar type stars, and, pulsar compan¬ 
ions asides, they have only been found orbiting such stars. 

* Partly based on observations made at the Observatoire de 
Haute Provence (CNRS) 

** Partly based on observations obtained with the swiss 1.2m 
telescope at the European Southern Observatory 
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This reflects to some extent an understandable desire to 
identify close analogs to our own solar system, which could 
perhaps contain life sustaining planets. Also, the selection 
function of the radial velocity planet searches has a rela¬ 
tively sharp optimum around spectral class G. Essentially 
all stars hotter than approximately F5 have fast rotation 
(Wolf et ah, 1982), so that it is impossible to measure their 
radial velocity to the ~10m.s“^ accuracy needed to detect 
planets. At the other end of the mass spectrum, most M 
dwarfs have slow rotation (Delfosse et ah, 1998a) and their 
velocity can be measured accurately, as we discuss below. 
Their luminosities however are much lower than those of 
solar type stars. At a given distance a much longer integra¬ 
tion time is thus needed to obtain a given radial velocity 
precision on an M dwarf than on a G dwarf. All planet 
search programs have thus understandably concentrated 
on solar type stars. 

G dwarfs however only represent a small fraction of the 
disk stellar population, with the lower mass M dwarfs out¬ 
numbering them by about an order of magnitude (Gliese 
& Jahreiss, 1991). It is thus likely that most planets in 
our galaxy orbit stars whose mass and luminosity are sig¬ 
nificantly lower than the Sun’s (Boss, 1995), unless some 
as yet unidentified physical process restricts planet forma¬ 
tion to the environment of sufficiently massive stars. It is 
clearly important to establish whether such a mechanism 
exists. 

For the last three years, we have been monitoring the 
radial velocities of a nominally volume limited sample of 
125 nearby M dwarfs. The two main goals of this large 
observing program (~30 nights/year) are to establish the 
controversial (e.g. Kroupa, 1995, and Reid & Gizis, 1997, 
for two contrasted views) multiplicity statistics of field M 
dwarf systems, and to pin down the still uncertain mass- 
luminosity relation at the bottom of the main sequence. 
Delfosse et al (1998b) present preliminary results for the 
stellar companion search, with 12 new components found 
in these nearby M dwarf systems, including the third de¬ 
tached M dwarf eclipsing binary (Delfosse et ah, 1998c). 
A byproduct of this program, related to the angular mo- 
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mentum dissipation of very low-mass stars, is described in 
Delfosse et al (1998a). 

Even though this was not the main focus of the pro¬ 
gram, we also realised from the start that for most of 
these stars we obtain radial velocity precisions which are 
sufficient to detect giant planets, if any exists around 
them. We present in this letter the first such detec¬ 
tion, around G1 876 (BD-15°6290, LHS 530, Ross 780, 
HIP 113020), a V=10.2 M4 dwarf (Reid et ah, 1995) at 
d = 4.702±0.046 pc (ESA, 1997). 

Delfosse et al. (1998a) present in detail the observed 
sample, while Delfosse et al. (1998b) discuss the observ¬ 
ing and analysis technique at length. We therefore only 
briefly summarize this information in section 2. We then 
proceed to discuss in section 3 the radial velocity detec¬ 
tion of the planetary companion of GL 876. In section 4 
we consider the implications of this detection and suggest 
some possible follow-up observations. 

2. Observing program 

The sample contains the 127 M dwarfs listed in the third 
edition of the nearby star catalog (CNS3 preliminary ver¬ 
sion, Gliese & Jahreiss, 1991) with a distance closer than 
9 pc, a B1950.0 declination above -16 degrees, brighter 
than V=15, and without a close much brighter primary. 
Observations have been carried out since September 1995 
with the ELODIE fiber-fed spectrograph (Baranne et ah, 
1996) on the 1.93m telescope at Observatoire de Haute 
Provence (OHP). The R=42000 spectra are wavelength 
calibrated through simultaneous observations of a tho¬ 
rium lamp. Since June 1998 some southern stars have also 
been observed with the nearly identical CORALIE spec¬ 
trograph on the recently commissioned swiss 1.2m tele¬ 
scope at la Silla (Chile). CORALIE mostly differs from 
the older ELODIE instrument by its spectral resolution 
of R=50000, better sampling of the spectrograph PSF by 
the CCD camera pixels, and significantly better tempera¬ 
ture control. The first indications are that these modifica¬ 
tions together result in a substantially improved intrinsic 
stability (Queloz et ah, in preparation). 

The extracted M dwarf spectra are analysed through 
cross-correlation with a binary (O/I) template constructed 
from an observed ELODIE spectrum of Barnard’s star, 
G1699 (Delfosse et ah, 1998b). For slowly rotating stars 
the resulting velocities have internal standard errors (pho¬ 
ton noise plus low level uncalibrated instrumental insta¬ 
bilities) which typically range from 10-15 m.s“^ for bright 
M dwarfs (V ^ 10) to ~75 m.s“^ at the magnitude limit 
of the sample. For Cl 876 (V=10.2) typical standard er¬ 
rors are 10 to 20 m.s“^, depending on airmass and see¬ 
ing conditions. Magnetic activity is common in M type 
dwarfs, and may further degrade the measurement accu¬ 
racy (Saar et ah, 1998). This potential error source is still 
uncompletely characterised for M dwarfs, but for slowly 
rotating stars (V sini < 3km.s~^) we can already bound 


Table 1. Orbital elements of GL 876. 


Element 

Value 

St. Err. 

Unit 

P 

60.97 

.19 

Days 

T 

2450661.7 

1.5 

Julian Days 

e 

.336 

0.019 


w 

5.2 

4.8 

deg 

Ki 

248.0 

6.6 

m.s"^ 

Vo 

-1.902 

0.006 

km.s"^ 

ai sin i 

0.00131 


AU 

f(m) 

7.810"® 


Mq 

O-C(CORALIE) 

23 

(rms) 

m.s"^ 

O-C(ELODIE) 

16 

(rms) 

m.s"^ 


it to (Tvr ^ 20m.s“^ for our cross-correlation analysis with 
the M4 binary template. Within the brighter two thirds of 
the sample, a conservative assumption at the present time 
is that we will detect all variables with semi-amplitudes 
larger than 40m.s“^. Assuming for illustration a 5 years 
period, this corresponds to a 1 Jupiter mass (Mj) planet 
orbiting a 0.25 M© M4V primary (at 1.8 AU), or to a 
2 Mj planet orbiting a 0.6 M© MOV primary (at 2.5 AU). 

3. A planet around G1 876 

Since planet detection was not initially emphasized in the 
observing program, its sampling strategy is not optimal for 
detection of low amplitude variations on timescales shorter 
than a few years. Cl 876 was observed once at each ob¬ 
serving seasons in 1995 and 1996 and its velocity variations 
became apparent from the three observations obtained in 
late 1997. It was then marked in our program lists as a 
variable. This low declination source however became un¬ 
observable from OHP before we could gather more data 
and determine its orbit. The commissioning of the swiss 
1.2m telescope at la Silla and its CORALIE spectrograph 
in June 1998 provided the first opportunity to obtain 3 
additional measurements of this southerly source, which 
allowed to finally determine its orbit. These observations 
were obtained within two weeks of the first light of this 
telescope, providing an encouraging indication on its po¬ 
tential for planet discovery. An end of night measurement 
from OHP at a large airmass provided a confirmation on 
June 22, just in time to confidently announce the discovery 
at the lAU “Precise Stellar Radial Velocities” conference 
(Victoria, Canada, June 21'** to 26*^). At this conference 
we learned from G. Marcy that his group independently 
identified the orbit of Cl 876, with orbital elements com¬ 
patible with our own determination. Weather permitting, 
we have since then attempted to observe Cl 876 at most 
every three nights, and often every night. The 1998 data 
therefore dominate the orbital solution. 

The orbital solution is given in Table 1. Preliminary 
solutions included a velocity zero point offset between the 
northern (ELODIE) and southern (CORALIE) datasets 
as a free parameter. The two systems were found to be 
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Fig. 1. Combined ELODIE and CORALIE radial veloci¬ 
ties for G1 876. The solid line is the radial velocity curve 
for the orbital solution. 

entirely consistent, and this parameter was thus held fixed 
to zero for the final solution. Figure 1 shows the individ¬ 
ual radial velocity measurements as a function of orbital 
phase (the 16 orbital periods elapsed since the first mea¬ 
surement make unpractical a display as a function of time; 
we however have essentially continuous coverage of one pe¬ 
riod in June and July 1998, excluding any possible spectral 
alias). The orbital period is two months and the velocity 
semi-amplitude is ~250 m.s“^, over 10 times the stan¬ 
dard error of one radial velocity measurement. The radial 
velocity curve implies a moderate but highly significant 
eccentricity of e=0.34. 

The large amplitude and moderate period of the ra¬ 
dial velocity variation argue strongly for orbital motion as 
its cause. An integration of the radial velocity curve im¬ 
plies a minimum physical motion of ^0.5 Rq. This varia¬ 
tion is about twice the radius of an M4 dwarf (Baraffe & 
Chabrier, 1996; Chabrier & Baraffe, 1997), excluding pul¬ 
sation as a possible explanation. G1 876 in addition only 
has low level photometric variability, and indeed happens 
to be one of the standard stars of the original UBV sys¬ 
tem (Johnson & Harris, 1954). It is actually variable, but 
with low rms amplitudes of 13 mmag at V, 9 mmag at 
R and 6 mmag at I (Weis, 1994). The photometric vari¬ 
ations don’t phase well at 61 days and appear consistent 
with a BY Dra type variability (Sasselov and Gody, pri¬ 
vate communication). Since G1 876 is a very slow rotator 
(Vsinz <2km.s“^, Delfosse et ah, 1998a), rotational mod¬ 
ulation of such low level surface inhomogeneities cannot 
explain its large radial velocity variations. Densely sam¬ 
pled photometry would on the other hand be of consider¬ 
able interest to establish the stellar rotation period. 

The mass of G1 876 unfortunately contributes some 
uncertainty to the minimum mass of its companion. As a 
consequence of H 2 recombination in the photosphere and 
the deepening convection for lower mass stars (Kroupa et 
al. 1990), the luminosity does not drop nearly as quickly 
per unit mass for mid-M dwarfs as it does for both higher 
and lower mass stars (Henry & Me Garthy, 1993), and it 


has a stronger metallicity dependence. Between ~0.50 Mq 
and ~O.18M0, Mass-Luminosity relations therefore have 
both shallow slopes and large intrinsic dispersions (Henry 
& Me Garthy, 1993). G1 876 thus belongs to a spectral 
type range where the mass of a single star is poorly 
constrained by its observable characteristics. Taking at 
face value either the solar neighborhood observational 
mass-luminosity relation of Henry & Me Garthy (1993) 
or the solar metallicity models of Baraffe et al. (1998), 
the absolute magnitudes of G1 876 (My = 11.81, Mj=7.56, 
Mz/=6.96, Mk=6.70, Leggett (1992) and ESA (1997)) im¬ 
ply a mass of 0.30±0.05 Mq for G1 876. As an illustration 
of possible uncertainties however, Delfosse et al (1998c) 
measure M = 0.432 ±0.001 M 0 and My = 11.7±0.2 
(My=11.81 for G1 876) for the brighter star in GJ 2069A, 
an M3.5V eclipsing binary which is probably super-metal- 
rich ([M/H]~±0.5). From its position in colour-colour di¬ 
agrams (Leggett, 1992), and from the relative depth of 
its cross-correlation dips with several binary templates 
(Delfosse et al., in preparation), G1 876 is probably more 
metallic than the sun, though not as much as GJ 2069A. 
We adopt a mass of 0.3 Mq for the rest of the discussion 
but warn that it is uncertain by perhaps 30%. The min¬ 
imum semi-major axis and planetary mass which result 
from the orbital solution are then asini = 0.20 AU and 
M 2 sinz = 2.0Mj. They respectively scale as Ms and Ms. 

4. Discussion 

The orbital elements of Gliese 876b are worth noting. In 
spite of the low mass of G1 876, the ice-condensation radius 
at the time of planet formation around this star is ~4 AU, 
only 20% lower than around a solar type star (Boss 1995). 
Once again the orbital semi-major axis (a = 0.20 AU) is 
thus much smaller than the expected minimum radius for 
giant planet formation, and some orbital migration must 
have occcured. However the observed orbital separation is 
also 4 times larger than the measured semi-major axes of 
51 Peg, r Boo and v And (0.04-0.05 AU). The excess of 
planets with such small semi-major axes is believed to re¬ 
sult from outward torques which counteract at short dis¬ 
tances the inward torque induced by the interaction of 
the planet and the protoplanetary disk (Lin et al. 1996, 
Trilling et al. 1998). These torques only become effective 
at separations significantly smaller than 0.20 AU, and can¬ 
not have played a signficant role for the G1 876 system. 
It is also interesting to note that the orbit of GL 876b 
is eccentric (e=0.35), while interaction with an accretion 
disk is expected to damp any significant orbital eccen¬ 
tricity of a planet (Goldreich and Tremaine, 1980). Sev¬ 
eral mechanisms may be called in to explain the large or¬ 
bital eccentricities of giant extrasolar planets, but in the 
present case the most interesting possibility is probably 
the chaotic interaction of several giant planets (Weiden- 
schilling & Marzari, 1996; Rasio & Ford, 1996; Lin & Ida, 
1997). A frequent final result of such a strong gravitational 
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interaction is a planetary system with a single giant planet 
at a moderate semi-major axis, in an eccentric orbit. This 
could thus simultaneously explain the semi-major axis and 
the eccentricity. 

Contrary to all previously confirmed planets around 
main sequence stars, G1 876b orbits a star which is very 
different from our Sun, showing that planetary systems 
form around stars of widely different types. G1 876 is 
much less massive than the Sun, ^O.SMq, and at most 
only ~150 times more massive than its planet. Its ra¬ 
dius is only three times as large: the radius of G1 876 is 
~0.3 i? 0 (Chabrier & Baraffe, 1997), while that of Jupiter 
is 0.10 solar radii. G1 876 is also much cooler than the Sun, 
and much less luminous. From the observed I-K and V-I 
colours (Leggett, 1992) its effective temperature is 3100 
to 3250 K (Leggett et ah, 1996), compared with 6000 K 
for the solar surface. From its absolute V magnitude and 
the bolometric correction of Delfosse et al. (1998a) it bolo- 
metric magnitude is 9.42, corresponding to 1.35 10“^ L©. 
Even though the planet of G1 876 is twice closer to its star 
than Mercury is to the Sun, the stellar flux at G1 876b is 
only ten times the solar flux at Jupiter, and lower than 
the solar flux at Mars. The apropriate albedo for G1 876b 
is unclear, and, by analogy with Jupiter (e.g. Podolak et 
al., 1993), its thermal balance of G1 876b may also include 
a substantial contribution from an internal heat source. A 
detailed evaluation of its effective temperature is thus be¬ 
yond the scope of the present letter, but G1 876b is clearly 
much too cold to possibly sustain liquid water above the 
1 bar level. 

G1 876 is closer to us than all other stars orbited by 
known extra-solar planets, by at least a factor of 3. At 
d=4.7 pc, G1 876 is the 40*^ closest stellar system to our 
Sun, and the 53’’'^ closest star. Since M dwarfs make up 
~80% of the solar neighbourhood population (Gliese & 
Jahreiss, 1991), it is only natural that the first member of 
this numerous class found orbited by a planet is a very 
nearby one, unless planetary formation would have se¬ 
lected against low mass stars. This discovery weakens such 
an hypothesis but improved statistics would obviously be 
needed for a reliable conclusion. 

This detection represents an opportunity to confirm 
a radial velocity detected planet through astrometry and 
determine its actual mass, or at the very least to set a 
lower limit which is firmly planetary. G1 876 is both at 
least 3 times closer to us than any other star with a de¬ 
tected planetary companion, and about 4 times less mas¬ 
sive (only 0.3 Mq instead of about ~ IMq for all previ¬ 
ous detections). Despite its relatively short orbital period 
of 61 days, the astrometric reflex motion induced by its 
~ 2 Mj companion is therefore unusually large by extraso¬ 
lar planet standards, with a minimum semi-major axis 
of 0.27 milliarcsecond for an edge-on orbit, and larger 
by l/sin(i) for more face-on geometries. The best single 
measurement precision of an astrometric observations is 
at present 1 milliarcsecond, with the FGS instrument on 


HST. A detection is clearly an ambitious measurement at 
this time if the orbit is seen edge-on, and it would need 
a very determined effort. A lower limit on the inclination 
of I sini| > 0.25 on the other hand will be easily obtained, 
and would already imply M 2 <8 Mj. In addition, these 
observations can be accomplished over the short timescale 
of one orbital period, only 2 months. 

Finally, it is interesting to note that the measurements 
of G1 876 obtained with the new GORALIE spectrome¬ 
ter on the 1.2 m telescope have residual 0-Cs as small 
as 22 m/s, for a V magnitude of 10.2. This discovery 
of a giant planet around a rather faint M4 dwarf illus¬ 
trates that this small telescope will powerfully contribute 
to the search for extrasolar planets. The application of the 
cross-correlation technique to the full wavelength domain 
(300 nm) of GORALIE compensates to some extent the 
disadvantage of the small telescope aperture (Baranne et 
al. 1996), and the nearly full-time availability of the tele¬ 
scope for planet searches will make future period identifi¬ 
cations much easier. 
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